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ABSTRACT 

We discuss properties of interpolating geometries in three dimensional gravity in the presence of 
a chiral anomaly. This anomaly, which introduces an unbalance between left and right central 
charges, is protected under RG flows. For this simple reason it is impossible to gap a system with 
such an anomaly. Our goal is to discuss how holography captures this basic and robust feature. 
We demonstrate the absence of a mass gap by analysing the linearized spectrum and holographic 
entanglement entropy of these backgrounds in the context of AdS 3 /CFT 2 . 
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1 Introduction 


One ambitious question that has long confounded field theorists is the classification of possible 
IR behaviour for a given UV fixed point subjected to arbitrary relevant perturbations. A, C 
and F theorems, whose various variations proved in 2d, 3d and 4d constitutes a powerful set of 
constraints on possible renormalization group (RG) flows in Lorentz invariant and unitary theories. 
The most notable success being Zamolodchikov’s theorem [1], and the recent proof of Komargodski- 
Schwimmer [2, 3]. 

The AdS/CFT correspondence has brought many new insights in understanding this classifi¬ 
cation. One radical feature in the correspondence is that different characteristics and phenomena 
in a quantum field theory—such as RG flows or entanglement entropy between different regions— 
manifest themselves geometrically in the dual gravity description. Constraints on RG flows have 
been studied exentsively holographically^ and insights learned in this process, such as the con¬ 
nection to entanglement entropy [9], has been crucial in understanding the above theorems. In 
particular, its influence in the three dimensional F-theorem [10, 11, 12] is without doubt. 

There are other features that constrain RG flows in quantum field theory. Most notably, quan¬ 
tum anomalies play a significant role. For example, it was argued by ’t Hooft that ABJ type 
anomalies of continuous global symmetry groups are necessarily conserved along a RG flow, i.e. the 
UV and IR theory suffer from the same amount of anomaly [13]. Very generally, quantum anoma¬ 
lies have shown up as the most important defining signature of the boundary states of different 
phases of matter found in condensed matter, most notably the classes of state named “symmetry 
protected topological” (SPT) phase [14, 15, 16]. In that context, it is also argued from very general 
grounds that boundaries of these exotic phases are robust against developing mass gap via adding 
any symmetry preserving perturbations in the infrared [14, 17, 18]. 

Even though it is understood how to design a holographic system that can flow to a gapped 
system in the IR, it has not been established how quantum anomalies are robust and protected 
geometrically. It is thus a curiosity how such “un-gappability” presents itself geometrically, where 
one might expect some restrictions over admissible background solutions. Our aim is to capture the 
“un-gappability” holographically and quantify how geometry knows about the powerful theorems 
for quantum anomalies. 

The current paper is a step toward extracting a geometric signature of quantum anomalies. 
While it would be interesting to study similar effects in higher dimensions, the discussion here will 
be limited to two dimensional QFTs and their three dimensional holographic duals. 

^There is an overwhelming amount of literature on this subject. We refer the reader to [4, 5, 6] for early work on 
the subject; see [7, 8] and references within for a more current status on the subject. 
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1.1 Chiral anomalies 


The simplest and most elegant example in field theory where a protected anomaly arises is the 
chiral anomaly in 1+1 dimensions. This is a theory where the OPE of the stress tensor is given by 


T{w)T{0) ~ 
f{w)f{0) ~ 


cl/2 T{0) dT{0) 

d ' 9 ' ' 

cr/2 f ( 0 ) dfjO) 

—d ' —9 ' — ' 


( 1 . 1 ) 


with T{w) = 27rTwu, and T(w) = 27 rr^^; moreover, its chiral nature comes from having c = 
Cl — cr ^ Q. Cl and cr are commonly denoted the left and right central charges. Whereas 
c = {cl + cr)/2 can flow according to Zamolodchikov’s theorem, it is also well known that c 7 ^ 0 
is protected.^ And for this simple reason, adding any symmetry preserving perturbation will not 
generate a mass gap in the infrared. 


The chiral anomaly provides a robust constraint on RG flows. But, how does the 2+1 holographie 
dual theory know about this constraint? Examples of gravity duals of chiral CFTs are already known, 
such as, most notably, topologically massive gravity theory in three dimensions (TMG) [21, 22, 23]. 
This is thus an ideal playground to observe how this protected anomaly plays out in the bulk 
space-time. 


The relevant comparison in holography is to contrast results of Einstein gravity coupled to 
matter versus TMG coupled to matter. In Einstein’s theory it is possible to craft geometries whose 
UV structure is conformal and the IR behaviour can be designed to be gapped. The surprise is 
that these backgrounds persist to be exact classical solutions in TMG! It must then be the case 
that, even though the background is unchanged, observables change dramatically after including a 
chiral anomaly in the bulk three dimensional theory. 

Given the expectation that the boundary is robust against arbitrary relevant perturbations, we 
will study perturbations around a generic RG geometry preserving boundary Lorentz invariance 
sourced by a scalar field with very general interaction potentials V{4>), and demonstrate how the 
boundary GET stays gapless all the way. This is in sharp contrast to the case of Einstein the¬ 
ory where one readily finds holographic RG flows leading to a gapped infrared theory. We will 
also confirm such expectation by studying holographic entanglement entropy. The standard Ryu- 
Takayanagi prescription [24, 25] is modified by this anomaly, which was developed in [26], and its 
modification is crucial to account for the expectation from the dual boundary theory. This is a 
non-trivial test of the proposal in [26], since the success shown here does not rely on symmetries 

simple argument goes as follows: modular invariance, in particular T transformations, imposes that 24c £ Z, 
and hence c cannot continuously vary. An argument that is more general goes along the lines of the analysis presented 
in section C for conserved currents. See also [19, 20] for a complete derivation. 
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of the background: it is a dynamical feature of the probe used to evaluate entanglement entropy. 
The effects of gravitational anomalies on entanglement entropy, on both sides of the holographic 
duality, have been recently generalized to higher dimensions in [27, 28, 29, 30]. 

Another approach would be to build holographic c-functions where the energy scale is parametrized 
by the radial direction. This was the approach used in [20], and they showed that their was no 
radial dependence on c. In our opinion, this approach carries some ambiguities since there is no 
robust definition of such c-functions (or at least the precise interpretation of it). Our derivations 
and conclusions do not rely on these definitions, nevertheless we are in perfect agreement with the 
results reported in [20]. 

1.2 Outline 

This paper is organized as follows. In section 2 we will review the criteria that allows for a given 
bulk solution to be interpreted as a RG flow that interpolates between a conformal UV fixed point 
and a gapped theory in the IR. The two key observables that we will use to signal a gapped 
IR are linearized spectrum and holographic entanglement entropy. The discussion in section 2 is 
valid for any dimension, but in later sections we will limit to three dimensional bulk theories. In 
section 3 we introduce a chiral anomaly by including a gravitational Chern-Simons term in three 
dimensions. For generic non-AdS backgrounds, we will show that the IR physics are dominated 
by the gravitational anomaly. This is clearly reflected by studying both linearized perturbations 
and holographic entanglement entropy. In the first two appendices we cover content related to 
conventions and details on the results related to section 3. And in appendix C we discuss how 
similar results holds for holographic theories that have unbalanced gauge anomalies. In this case 
the effects of a quantum anomaly on the RG flow is not reflected classically in the bulk geometry; 
instead we find that quantum effects in the gravitational dual have to be included to appropriately 
capture the “un-gappability” of the system. 


2 Gapped holographic systems: entanglement entropy and spec¬ 
tral analysis 


The class of holographic systems we will study are vacuum (i.e. zero energy, zero entropy) solu¬ 
tions which in the UV are conformal and the interior IR geometry has a non-trivial radial profile. 
Depending on this radial profile, one could infer that the IR geometry is dual of a QFT with a 
gapped spectrum. There are two common holographic routes to determine if the IR geometry 
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is gapped or not.^ The first route is to explicitly quantify the spectrum by studying linearized 
perturbations around the given background of interest, see for example [6, 33, 34]. The second 
route is to study the behaviour of the entanglement entropy, as in [35, 36, 34]. These are the two 
observables (spectrum and entanglement) that we will discuss throughout this work. In this section 
we will summarize the main features of these observables when the dynamics of the holographic 
background is governed by an Einstein-Hilbert action coupled to matter. 

Our starting point is to consider geometries which we parametrise as 


ds^ = 


rjijdx^dx^ dz^ 


+ 


z‘^f{z) 


i = 0,... ,d 


( 2 . 1 ) 


As z —)• 0, i.e. in the UV limit, the spacetime asymptotes to AdS^+i which requires f{z) —)■ 1. In 
the infrared limit, z —)> oo, we will take f(z) ~ z^ for some power n > 0. This geometry can be a 
solution of an Einstein-Scalar system 


i[9, 0] = {R - 8(v</.)2 - vi<p)) . 


( 2 . 2 ) 


Provided f{z) ~ z^, it is sufficient to have 


V{(j)) ~ exp(—4(I(/)) , 


(2.3) 


in which case the power n can be expressed in terms of 5 as 

n = id-l)6^. (2.4) 

The null energy condition requires n < 2d. In the following we will determine for which ranges of 
n the infrared geometry can be interpreted as dual to a gapped system. 

Before proceeding it is important to note that these solutions are singular in the interior: as 
z ^ oo there is a curvature singularity in (2.1) and the scalar field diverges. This is an apealing 
feature since it makes evident that the background is not a black hole (i.e. a thermal state), and 
hence it can be interpreted as the RG interpolation between two theories. For our purposes the 
singular nature of the geometry does not spoil the analysis: what is important is that we can 
properly quantify the observables which will be evident below. Moreover, for 0 < n < 2d, it has 
been argued extensively that the singularity can be resolved either by embedding into string theory 
or other methods [33, 37]. 

The simplest way to quantify the spectrum, is to consider a massless probe scalar 'I'(x^) on the 

^There are of course more ways to test if the IR point is gapped. For instance, the thermodynamics response of 
the system; see [31, 32]. 
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background (2.1). The resulting Schrodinger equation is 


si'll; + Ves{z{u))i; = k^ij; , V'eff(z) 



(2.5) 


where 

vl/(z, xn = , duiP = yW)dz'>P , (2.6) 

and k? = 'qijk’^k^; the radial coordinate u is defined appendix A. In the IR limit, we have f{z) ~ z"' 
and the effective potential becomes 

Feff ~ (2.7) 


For n > 2, as z —>■ oo, the potential I4ff diverges. Since I4ff also diverges in the UV limit, this 
suggests that we have an infinite well, and the spectrum is discrete in this well. Therefore, there 
should be a mass gap. For n < 2 however, vanishes as 2 ; —)• 00 , suggesting that it is a half well, 
and the spectrum is continuous and hence ungapped. The case n = 2 requires more care, and we 
refer to [33, 34] for a complete discussion. 

This simple example of a probe scalar gives us a range of n for when one should expect a gapped 
spectrum (2 < n < 2d) or not (0 < n < 2), and for many cases this probe is enough. However, it 
will be evident in later sections that probe fields are blind to the true nature of the system since 
they might not capture key aspects of the spectrum. A proper treatment requires studying the 
actual metric and background fields perturbations as done in e.g. [6, 33]. In appendix B.l we do 
this explicitly for a three dimensional version of (2.2). 




Figure 1: Schematic representation on the behaviour of minimal surfaces for a geometry that is not gapped 
(left picture) contrasted to a gapped background (right picture). The solid red lines correspond to connected 
surfaces, which are anchored at the boundary. The dotted red lines are disconnected surface, which dominates 

for R > Rmax- 

One could alternatively use holographic entanglement entropy as a probe of the spectrum: the 
topology of the entangling surface changes dramatically depending on the value of n, which is a 
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signal of a gap in the system [35, 36, 34], The argument can be summarised as follows. Consider 
the relation between the length of a strip at the boundary, which we denote R, as a function of the 
turning point of the minimal surface, denoted zq. For the background (2.1) this would give 


r^o 


R = 2 


dz 


^d-l 




( 2 . 8 ) 


For large zq we can replace f{z) ~ /o^"', and we will have 


R ~ 2zl + 0{zq ^)) , X = 




Jo Vl — u 

while for small zo, the leading effect is given by /(z) ~ 1, and we get 


2 ’ 


(2.9) 


R^2zo + O{zo) 


( 2 . 10 ) 


What these manipulations clearly illustrate is that if n > 2 then i? —)• 0 both as zq increases and 
decreases. This indicates that connected geodesics have a maximal value R'^^^ when n > 2, and 
hence for boundary intervals with R > R^^ the disconnected solution would take over. This is 
illustrated in figure 1. The absence of R in 5 ee for disconnected geodesics signals that the system 
is gapped in the IR. The analysis and bounds on n using entanglement is in complete agreement 
with that of the spectral analysis reviewed above. ^ 


3 Chiral anomalies in AdS3/CFT2 


A simple way to prevent a gap in the IR is to introduce a protected anomaly. In a two dimensional 
CFT this can be done by adding a chiral anomaly, i.e. by having an unbalance in left versus right 
central charge, cl / cr. Whereas the total central charge c = (cl + cji)j2 is a monotonic function 
in the space of RG flows, the difference c = cl — cr is protected: given a value of c in the UV, it 
remains intact in the IR. In this section we will show how three dimensional gravity protects c. 

Holographically it is known how to induce a non vanishing value of c in three dimensional 
gravity. Such systems are described by the action principle (2.2) with the addition of a gravitational 

^For d = 2 and n = 1 there is a small caveat to this line of reasoning. There are situations where the connected 
solution is always dominant, but as it —^ oo the length of the curve does not depend on R. This is simply due to 
the fact that the effective central charge asymptotes to zero in the IR, and it could be interpreted as a signal of a 
“marginal gap” in the spectrum. In this context marginal refers to the deformation inducing a very rapid power law 
decay of observables with distance, rather than exponential decay (which we would denote as a “hard” gap). 
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Chern-Simons term 


-^GCS 


1 


327rG3^ 


(rL 8 ,r"„ + \TiXgrK, 


(3.1) 


Gravitational theories that include both the Einstein-Hilbert term and /gcs are known as topologi¬ 
cally massive gravity (TMG) [21, 22, 23]. The role of this term is to induce a gravitational anomaly. 
In particular for a space-time that is asymptotically AdSs the gravitational Ghern-Simons term in¬ 
duces an unbalance in the central charges given by [38, 39, 40, 41] 


CL 


3£ 


(1 + 



CR 




(3.2) 


Our goal is to show that for finite /r it is impossible to have a background solution that could be 
interpreted in the IR as a gapped system. 

One might think that the obstruction in gapping the IR geometry would be to argue that a 
“gapped” background ceases to exist after including /gcs- This would have been a clean and 
elegant resolution, however, it is not the case. The contribution of /gcs to the equations of motion 
is given by 


Rfiv - = 8d^4>d^(j) - 45 ^^ 5 "(/>( 9 o (/> - {^) . ( 3 . 3 ) 

where is the Gotton tensor: 

C'm. = - \Rgfiu) = lie/’^VpRau + cr^pRap) . (3.4) 


Despite the complexity of this equation, it was noted in [20] that for all metrics of the form 


ds^ 


dz^ T]ijdx^dx^ 


+ 


z‘^f{z) 

, 'o + gijdx^dx^) , 
z{uY ^ 


(3.5) 


the Gotton tensor vanishes for arbitrary z{u). This is simply due to the fact that the metric (3.5) is 
conformally flat and thus C^u = 0. This means that if (3.5) is a solution of a two derivative action, 
it would automatically be a solution after adding /gcs- All backgrounds discussed in section 2 
fall into this category, and hence it remains a solution even when cl Y cr- Moreover, using the 
prescription of e.g. [20, 7, 8], the holographic c-function is given by 


2G3/(z)V2 


(3.6) 









which is unaffected by the gravitational Chern-Simons term since it is built out of traces of the 
equations of motion. This quantity suggests strongly the total central charge c = cl + c/j is 
decreasing monotonically along the radial ‘RG’ direction. But this is of course not enough to infer 
what happens to cl{z) and cr{z) independently; moreover the dual interpretation of (3.6) far from 
the fixed points is somewhat ambiguous. Still this raises some concerns; if we naively think that 
the results of section 2 still hold we would interpret (3.5) as a gapped system in the IR for a theory 
with Cl 7 ^ cr. But this cannot be true! We just argued that it is impossible to gap a UV system 
with unbalanced left versus right moving excitations. Holography seems to say the opposite. In 
the following we will resolve this puzzle by quantifying the linear bulk modes and holographic 
entanglement entropy in the presence of the gravitational Chern-Simons term. 

3.1 Linearized spectrum 

Our first step towards resolving our puzzle is to look at the linearized perturbations of the matter 
fields that support the background (3.5). The difficult portion is encapsulated in the Cotton tensor 
in (3.3). Prior literature on linearising includes e.g. [42, 43, 41]. Within this literature, it 
is very well known that TMG has negative energy excitations, unless we are at the chiral point 
[44, 45, 46]. This feature will still be present here, however we will not highlight this explicitly 
below. Our aim is to quantify if the spectrum is discrete or continuos in the IR geometry. 

For the purpose of describing the linearized equations we will use a radial gauge 

9piu = , (t) = (/>(o) + 6(1) , huf, = 0 , (3.7) 

where {gjiJ, (f)(o)) define the background solution, and {h^y,6(j)) are the small perturbations around 
the background. Details of the computation are presented in appendix B. It is convenient to write 
the background solution (3.5) as 

ds'^ = gj^Jdx^dx" = (du^ -|- rnjdx^dx^) , = z{u). (3.8) 

and we will use null coordinates for the boundary directions; 

w = —t + X , w = t + x , d = dw , d = du, . (3.9) 


The background scalar (I)q satisfies (B.7)-(B.8). 
The differential equations for h^y and 6(p are 


1 


(^dldhww 


2d‘^dhww + dd'^huiu) 


+ d hu)i 


e^dudhww + e^dudhww - 8e'^4>{o) d6(p = 0 , 
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'WW 


288“^ hww + 8^hww + 8‘^8h, 


'WW 


'WW 


'WW 


— e 


•'^8u8hwu, + e^8u8h^u, + 8e^(^(o) 864> = 0 ,(3.10) 


and 



These equations correspond to the (i, j) components of (3.3), and they take into account constraints 
conditions which arise from the {u,i) components in (3.3). In the above we defined 4>o = 8u4>o and 
A = 8uA{u). There is also a constraint which is obtained by tracing (3.3) and the linearized 
equation of motion for (f>. These equations can be found in appendix B. 


It is important to note that all terms multiplying in the above equations are independent 
of the conformal factor ig gimply because the Cotton tensor is invariant under Weyl 

transformations and our choice of coordinates makes this explicit. This will be key in the following 
subsections. 


Finally, it is possible to manipulate the above equations such that {hww-, <5</>) essentially decouple 


from {hww, hww), and hence {hww,5<i>) should be taken as source terms in the equations for hww 
and hww- The decoupled equations governing hww and 5(j) are given in the appendix B.2. While 
this decoupling is conceptually important, the details of the equations are not crucial and we refer 
the reader to the appendix for details. 

3.1.1 IR limit: conformal gravity 

Our aim is to understand the behavior of the modes and 6(p in the IR geometry. Along the 
lines of the discussion surrounding (2.7), the behavior of the modes in the interior dictates if we 
have bound states or not. 

In the notation (2.1), the IR geometry is given by z —>■ oo and f{z) ~ with 0 < n < 4 (the 
upper bound for n is dictated by the null energy condition). In terms of the radial coordinate u in 
(3.8), the IR geometry is given by 


u ^ z —i-cso if0<n<2 

u ~ log( 2 ;) —>■ oo if n = 2 , 


u ^ uo + z tiQ if 2 < n < 4 


(3.12) 


where uq is a positive constant.® In addition we have = z which diverges in the interior. 

And according to (B.7), the background scalar behaves in the IR as e"^(;(>(o) ~ which decays 

as z —)• oo for the range of n we are considering. 


®And for completeness, the UV boundary is located at u ~ z —>■ 0 with f{z) ~ 1. 
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Zooming into the IR geometry, we see that dominant terms in equations (3.10)-(3.11) are 


2d‘^dhu,w + ddP'hyjw + = 0 

2dd‘^hww + d^hwu, + d'^dhijjw^ = 0 

^du{s^hu,w - = 0 (3.13) 

In this limit we kept (9, d) fixed and we discarded terms like e"^du since they are subleading relative 
to dy. These equations are those relevant for conformal gravity, i.e. a gravitational theory described 
solely by /gcs-^ The Cotton tensor is insensitive to and dominant in the IR: this is already 
indicating that the behavior of the modes is independent of n which should be contrasted with 
the discussion in section 2 and appendix B.l. This is a rather powerful observation that has no 
counterpart in the absence of the gravitational anomaly. 

Let us solve explicitly (3.13). It is important to note that hy,w contributes as a non-homogenous 
term for {hy,y,, h^w) since we can solve for hw^ using the equations in appendix B.2. The homoge¬ 
nous solutions to (3.13) are simply 


- ( d'^dh - 

^ V 

- (dldh,j,y, - 


^^iV2dd u 


^^iV2dd u 


(3.14) 


which holds regardless of the details of the properties of the radial profile. In particular, both 
solutions are regular in the interior so we conclude that these modes have a continuous spectrum; 
moreover the details of the IR geometry are irrelevant. In appendix B.2 we discuss the IR behavior 
of 5(j) and hyjw] in contrast to (3.14) these modes are sensitive to the details of the IR geometry 
and their spectral functions can be gapped. But since there is a sector of the theory that is always 
continous, which (3.14) makes evident, we conclude that the whole system is not gapped in the IR. 


3.2 Entanglement entropy 

We have just shown that, despite the geometric properties of the background (3.5), in TMG all 
vacuum solutions have a continuous spectrum in the IR region. In the following we would like to 
see how this is captured by entanglement entropy. This will prove how physical observables are 
starkly different even though the background geometry is unchanged. 

In the presence of a gravitational anomaly, the Ryu-Takayanagi prescription [24, 25] is modified. 
The modification due to the gravitational Chern-Simons term can be interpreted as the functional 
of a massive spinning particle [26]. Holographic entanglement entropy is given by the minimum of 

®For a recent discussion on conformal gravity and a complete set of references see [47, 48]. 
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the functional^ 

+ sn • Vn^ , (3.15) 

where m and s are the mass and spin of the test particle, and As shown in [26], the 

replica trick hxes these coefficients to be 


m = 


1 

^ ’ 


1 

AGsiJ, 


(3.16) 


The precession of the spinning particle is described by two normal vectors to the worldline of the 
particle: and h^. The minimization of (3.15) is subject to the constraints that n^, and 

are mutually orthonormal: 


5 , 


constraints 


= J dr [Ain • h + X 2 n ■ v + Ash • v + A 4 (n^ + 1) + A 5 (h^ — 1)] , 


(3.17) 


where Ai(r) are Lagrange multipliers. Note that is timelike, and is a spacelike vector. The 
resulting dynamics of the probe is given by 

V[mu^ + VpVs^^P] = = -se^''\x , (3.18) 

which is known as the Mathisson-Papapetrou-Dixon (MPD) equation, and the constraint 


= 1 . 


(3.19) 


We are interested in solving (3.18), and evaluating (3.15), for backgrounds of the form (3.5). 
The first simplification in this case comes about by noting that the system preserves translation 
symmetries along x* = {x,t). This implies that (3.18) contains two conserved quantities - the 
momentum along x and t- which are 


Pj = rne‘^\^+ se^^{v^fdr(^^^ , (3.20) 


with Px and Pt the constants of motion. It is also useful to note that (3.20), together with (3.19), 
implies 


sdr{e\^) = PxV^ - PtV^ . 


(3.21) 


'For any vector we define, 


dr dr 

Throughout this section, r is an affine parameter along the worldline of the curve, i.e. 


= 1 . 
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The remaining component of the MPD equation (3.18) is the radial direction which reads 

m [9,(eV) + - {v^f)duA] - s e-^dr {e^^{v^drV^ - v^drV^)) = 0 . (3.22) 

Note that this equation is not independent: by taking derivatives of (3.20) combined with (3.19) 
gives (3.22). 

3.2.1 No static solutions 

Before discussing the solutions to (3.20)-(3.22) that will be anchored at the boundary, it is in¬ 
structive to highlight a general property of the MPD equations. Consider taking a static ansatz, 

i.e. 

u* = ^ = 0 ^ (u“)2 + = e-2^ . (3.23) 

dr 

In this case the MPD equations (3.20) reduce to 

P^ = m , Pt = -s e^^'^^v^d^A{u) , (3.24) 

where we used (3.23). Assuming that both s 7 ^ 0 and m 7 ^ 0, it is clear that (3.24) is only consistent 
with both Px and Pt being constant if = u~^. This choice of A{u) corresponds to having 

exactly AdSa as a background and, moreover, would follow the usual geodesic path. The other 
obvious possibility is to take = 0 in (3.24) but this would correspond to a solution that is 
disconnected -the endpoints of a single curve are not anchored at the UV boundary. 

Hence, for a general background with 7 ^ u~^ and (m, s) non-zero, the choice u* = 0 will 
not allow connected solutions to the equations of motion of the spinning particle. One might worry 
about this issue: it has been argued extensively in the literature that for a static background the 
Ryu-Takayanagi curve must lie on a constant time slice [49, 50]. The MPD equations clearly violate 
this restriction. Even though the background (3.5) is static, the gravitational Chern-Simons term 
is not invariant under t —>• —t. We suspect that due to the lack of parity invariance of the theory 
there is no good reason to restrict the curve to lay on a space-like slice. 

Finally, we note that for = u~^ and u* = 0, the conserved quantities satisfy 


mPt = sPx ■ 


(3.25) 


3.2.2 Connected solutions: perturbative analysis 

We have established that connected solutions to the MPD equations must have u* 7 ^ 0 for non-AdS 
backgrounds, which makes the task of building connected solutions much harder. We will first build 


13 


these solutions by perturbing the MPD equations in s (the spin of the probe). Our goal in the 
following will be to hnd a relation between the zq (the deepest point of the path) and the size R of 
the boundary region; the aim is to show a breakdown of the relations ( 2 . 8 )-( 2 . 10 ). 

To solve perturbatively the MPD eqnations it is useful to use the radial coordinate as the proper 
length of the curve. In this gauge we have 

= v'^du , v{u) = = (1 - {dutf + . (3.26) 

Using this parametrization (3.20) becomes 

Px = m e^v duX + s v^d'^t , 

Pt = me^v dut + sv^d^x . (3.27) 

We are interested in solutions to the MPD equations that are anchored at the boundary, i.e. the 
boundary conditions we will impose are 

R = Ax = x{uf) — x{ui) , At = t{uf) — t{ui) = 0 , ttj, nj —)• 0 . (3.28) 

Taking s as a small parameter, we expand the curve as 

OO OO 

x{u) = xo(u) + ^ s^Xi{u) , t{u) = to(«) + ^ sHi{u) , (3.29) 

i=l i=\ 

and for the conserved charges 

OO 

Pt,.=p^+ ■ (3-30) 

i=l 

Using (3.29) and (3.30), we will solve the MPD equations (3.27) order by order in s. Note that the 
spin s is not quantized since the probe is a classical particle. From (3.16), small s and m fixed is 
equivalent to taking p large and G 3 fixed. 

We start at order the solntions are just those satisfying the geodesic equation. In accordance 
to the boundary conditions (3.28), at this order we will use the static geodesic: 

Po = 0 , to{u) = 0 , (3.31) 

and 

^ x| ’ ^0 = ■ (3-32) 

yj Zq — z{uY 

Looking ahead, we re-introduced z{u) = which brings the line element (3.5) to the form 
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(2.1), and zq which will be turn out to be the deepest point of the curve. In this notation duZ = 
f^^^{z) with f{z) defined in (2.1). 

At linear order in s, the first equation in (3.27) receives no correction due to the spin coupling 
-i.e. the last term in (3.27). For this reason we set 

pf = 0, duxi = 0 , (3.33) 


and hence at linear order in s we find no correction to R. The second equation in (3.27) gives at 
linear order 


duh 


— {p\duXo + duVo) 

Po 


z 

pW4 - 


fMH] 

Zo j 


(3.34) 


From here it is clear that the solution is not static: we cannot set ti to be a constant without 
drastic consequences. Imposing (3.28), which sets Ati = 0, gives 


Pi = 


Axn 


(3.35) 


where we used that ADq = vo{uf) — vo{ui) = 2; Axq is the length of the boundary interval in the 
geodesic limit. Note that (3.35), and in general t{u), can be easily obtained from (3.21). For our 
boundary conditions (3.28) the exact relation is 


Pt = 


sAv 


(3.36) 


Moving on to second order in s, the MPD equations give for the temporal component 

dut2 = ^duxo , (3.37) 

Po 

where we made use of (3.33). Since we require At = 0 this sets = 0) hence without loss of 
generality we set t 2 = 0. For the spatial component, the quadratic correction in s is 


m 

ziu) 


duX2 




1 


2pgDo 


-2.9J + 


til 

zo 


+ 


(^o)" 


^ ~ (Pil ) +P2 


1 

2Pb 


. (3.38) 


In the first line we already made use of (3.33), and from the first to the second equality we used 
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(3.34). The value of pf is 


= . P.38) 

which is determined by demanding convergence of X 2 when integrating (3.38). Replacing gives 


dzX2 = 


2rin?ZQVQ 

From here the quadratic correction to Ax is 
dq q 


zdzf{z) + {zqp\ - f{zf/‘^f + jJ^ifiz) - fizo))^ ^ 

[xor 


\ f{zy/^' 


(3.40) 


A ^0 

Ax2 = ^ 


Jo fizoqV^^ (1 - . 


qdgfizoq) + {zop\ - f{zoq)^^‘^f + 

(1 - q^) 


(3.41) 


where we defined q = zj zq. In the small zq limit, which means that we are near the AdS boundary, 
we have f{zoq) ~ 1 and the interval size is not corrected: Ax 2 ~ 0. This simply implies that near 
the boundary the curve is well approximated by a geodesic. 

In the opposite regime, i.e. large zq limit, we have f{zoq) ~ fo{zoq)'^ and the interval is corrected 


R = Axo + Ax 2S^ + O(s^) 

~ ^ ^ ^((/o^o) ^ (/5^o-\/7o^ + 0((/o2:o )^'^^)) + O(s^) ) (3.42) 

where for each order in s we are writing the leading correction in zq. These corrections are our first 
indication that there is a dramatic effect as we turn on the spin: the dependence of the interval 
size with zq is dramatically changed since now for n > 2 there is no bound on R when zq is large. 

As we go higher in perturbation theory it is simple to see that the structure of the corrections 
are even in s for x{u) and odd for t{u). Moreover, a simple analysis of the structure of these 
corrections will imply that as we explore the limit zq large we will find 

Ax2i ~ . (3.43) 

One might worry about the convergence of these corrections for arbitrarily large zq and s < 1, 
and this is legitimate concern. Our perturbative expansion in s breaks down in the IR geometry, 
since we are neglecting derivatives in 2 ; that become more important as zq is larger. The formulas 
above just illustrate that we can construct connected solutions, and that R receives a non-trivial 
correction. In the following we will investigate the IR behavior of the solutions. 

®q: and /3 are 0(1) numbers coming from the integral over q. These numbers are not important for the general 
scaling with zo and are sensitive to the approximation chosen for the integral so we will not specify them further. 
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3.2.3 Connected solutions: IR behavior 


The IR corrections we found in the regime s m come with generically large coefficients as it 
is reflected in (3.42). For the purpose of convergence of this series, we will discuss the opposite 
regime: s ^ m. 

Taking m = 0 reduces the MPD equations (3.27) to 

Px = sv^d^t , Pt = sv^d'^x . (3.44) 

Two remarkable things happen in this limit. First, the equations do not depend on the warping 
factor of the metric: Hence any solution we find in this regime is insensitive to details of the 

interior geometry. Second, in this regime me can solve the equations exactly. And related to this, 
in this limit we can recover static solutions: it is perfectly consistent to set t{u) = 0 which just sets 
Px = 0 which is not possible if both m and s are non-zero. 

On top of all of these elegant simplifications, there is a more interesting reason why to consider 
the limit s P> m. In section 3.1.1 we argued that the IR limit of the linearized equations suppressed 
terms like e^du versus d^. Looking at (3.27) this would imply that for curves that explore the deep 
infrared geometry, the spin coupling terms dominates over the mass coupling. The limiting case 
in (3.44) hence governs the IR dynamics, which as expected agrees that conformal gravity governs 
the long distance dynamics of the system. 

In order to find solutions in this regime, it is convenient to rearrange (3.44) as 

Pt dlt = Px dlx , s^v^dlv = ^ (3 45 ) 

The solutions to (3.45) are rather simple to obtain in general form. However, we are interested in 
solutions satisfying (3.28). This sets Px = ^ and dut = 0. For the spatial direction we find 

duX = — 7 =™=== , h{u) = Ptu + hi . (3.46) 

y — h{u)^ 

with hi a constant. The turning point of this curve is defined by h{uo) = s; uq is the analogous of 
zq used in the previous subsection. For this solution we find that 

7“° s h^ 

R = 2 duX = — = uo + — . (3.47) 

JO 

The width of the spatial interval is finite and obviously independent of the details of the geometry 
in the interior, i.e. there is no appearance of A{u). This should be compared with the geodesic 
solution using this parameterization: uq is defined by = m/Px (see (3.32)), which makes 
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the dependence on the warping factor explicit. 

The next step would be to match the solutions for small and large s. It is important to emphasize 
that the small s expansion assumes that the subleading corrections (3.43) are small which implies 
that Axo (or alternatively zq) cannot be arbitrarily larger than the the mass gap (controlled by / o ). 
For zo arbitrarly large the perturbation in s breaks down. Similarly, for large s we are neglecting 
terms in the MPD equation that can become dominant for short intervals. Still, there will be a 
matching region where Axq is of the order of the mass gap, and numerically one should be able 
find overlap of both regimes in s. As the matching is done, the integration constants (such as hi) 
in the IR will be fixed by requiring continuity with the UV. We will leave this for future work. 

3.2.4 On-shell action of the spinning probe 

We have shown, at least perturbatively in s, that we can design connected solutions to the MPD 
equations for arbitrary R. In the following we will quantify how these solutions modify the behaviour 
of entanglement entropy in the presence of the gravitational anomaly. 

We will evaluate (3.15) up to order for a connected solution. In order to do so, we need to 
build the normal vectors to our curve and impose appropriate boundary conditions. Recall that 
the normal vectors satisfy the following orthonormality constraints: 

5'constraints ^ j ^ + '^2n • V + Ash • V + A4(n^ + 1) + A 5 (h^ - 1)] , (3.48) 

and that the tangent vector satishes (3.19). 

It is rather simple to build normal vectors to = (u“,u*,u^) for the background (3.5). For 
instance, two mutually orthonormal vectors are® 

= Fvv\ Fvv^^ . (3.49) 

Here is Q is timelike and Q is spacelike; F~^ = ((v^)^ — (u*)^) which adjusts the normalization 

appropriately. The vectors {Q^,Q^) are however not very convenient to evaluate (3.15) since they 
don’t a priori satisfy any reasonable boundary conditions. Following [26], it is natural to require 
that the timelike normal vector should point along the CFT time direction at the endpoints of the 
curve. This allows us to compare the orientation of normal vector at each endpoint as we transport 

®Our convention for the epsilon tensor is 

! - txu 1 

stxu — y~9 1 e —- 1= ■ 
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them along the curve. 

To implement unambiguously any precise choice of boundary conditions, it is useful to construct 
normal vectors that are Fermi-Walker transported along the solution. More concretely, we want to 
first consider vectors that are orthonormal to and also satisfy 

= -{qaVv°)v^^ , = -{qaVv^)v^ . (3.50) 


The advantage of this parametrization is that we can easily evaluate (3.15) for our choice of bound¬ 
ary conditions. Using (g^,g^), it is convenient to write 


= cosh(7/(T))g^ -|- sinh(r/(r))g^ , , 

where 7?(t) is a function that is adjusted such that at the endpoints^^ 


(3.51) 


'Hi — — (^4)cFT ■ 


(3.52) 


Given (3.51) and (3.50), the anomalous portion of (3.15) reduces to 


5, 


= s 


dr h ■ Vn 


= s {r]{Tf) - riin)) 

V q{n) ■ Hi - q{Ti) ■ m 


(3.53) 


The appeal of this equation is that we actually don’t need to solve for ??(t). We just need to build 
{q^,q^) and evaluate the inner product with (3.52). 

To build (g^,g^) it is actually useful to use (3.49); in particular, it is rather clear that each 
basis of normal vectors should be related as follows 


g^ = cosh(/i(r))(5^-|-sinh(/i(r))Q^ , 

g^ = sinh(/i(r))(5^ -|- cosh(h(r))Q^ . (3.54) 

The advantage is that we reduced the undetermined number of variables, since we only need to 
solve for /i(r) such that (3.50) is satisfied. And this is a straight forward procedure: plugging (3.54) 
in (3.50), the resulting equation for /i(r) is 

d^h{T) = -v^v^) 

= T («”'-+ f) ■ 

^°Here we define (9t)cpx as a vector that points in the time direction at the boundary, however it is normalized to 
ensure that = —1 with respect to the bulk metric. 
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To obtain the first line, we just made use of (3.49); to obtain the second line we used the MPD 
equations (3.20). From here, we find that (3.53) becomes 

/ (y^ _ y yt\ „ \ 

Sanom = •§ log ( -j— -® / drOrhir) . (3.56) 


With these results we can now evaluate (3.15). Including all contributions we find 





= 2m 


I 


fii/u' u “T ^anom 

dz 1 , {{v^ — vv^)f\ dz 




/(Z)V2 


V 


(jjX _ ^ yfj. 


10 


/(Z)V2 


fTMG) 

At this stage, it is important to emphasize that we have made no approximations so far in . 

The expression above is exact and correctly implements the boundary conditions (3.52). Also note 
that in the first integral we explicitly include a UV cutoff e since the proper distance is divergent. 
However, the last integral in (3.57) does not carry divergences as z —)• 0; this will be evident shortly. 

We now proceed to evaluate (3.57) for our perturbative solution. The first integral, which 
measures the proper length, gives 


= 2m I" ' 


zfizy/"^ 


2 : ^zl - z2 


+ 


p dz z{zl{f{z)-f{zo)) + z{zl-zYd,f) , 


L 


Jo f{zy^‘^ 


Zo(zg - z2)3/2 


+ 0(^).(3.58) 




And for the terms in 5anom we obtain 


O r dz « dz 

Jo f{zy/^^ mi /(z)V 2 ' 


+ 0(4) 


m'^ 


(3.59) 


and 


s log 


(u* — VV^)f 


(^yX _ ^ yfj. 

Combining these three contributions gives 


2zo 


2s2 

m V Axo 


- 1 


= 2m 


r 


dz zq 1 2s^ / 2zo 

/(-z)02 z ^z^ - z2 m VAxf 


- 1 


f zzoifjz) - /(zq)) 
Jo /(z)V2 I (^;2 _ ^2)3/2 


7^° dz 

+— 


+ dzf- 


2z^ — Zq 


Zo 


^4- 


rri'^ 


(3.60) 


+ 0( ), (3.61) 


fTMG ^ 

In the limit zq <C 1 where /(z) ~ 1, one immediately finds 5^^ reduces to the geodesic 
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length. This is of course expected as we are only probing the geometry near the boundary. The 
other interesting regime is for large zq where we use f{z) ~ foz'^. We find 

+ , (3.62) 

where “...” represents corrections both in and in I/zq- We immediately see that the order 

drastically changes the behaviour of the entanglement entropy. More importantly, by trading zq 
fTlVTC ^ 

with R via (3.42), S'gg is always dependent on the size of the interval. This is generic for gapless 
systems as expected. 


4 Discussion 

In this work we discussed properties of RG backgrounds in three dimensional gravity in the presence 
of a gravitational anomaly. Our main result was to show that the infrared geometry could never 
be interpreted as the dual of a gapped system, as expected from general QFT arguments. The 
geometry did not make evident that a gravitational anomaly is protected; the analysis presented 
here is an attempt to make these universal features robust and in accordance with the RG theorems 
in the dual theory. 

The most straightforward computation that captures how a chiral anomaly is protected is 
the spectrum of fluctuations. We showed that, regardless of the radial profile of the background 
geometry, the linearized spectrum of gravitational perturbations has a continuous sector deep in 
the IR. This can be a rather cumbersome task; a simpler holographic observable is entanglement 
entropy which we discussed at length. 

As expected, holographic entanglement entropy behaved in accordance with RG theorems: we 
showed that connected solutions are always present in any vacuum solution in the presence of the 
gravitational anomaly, which makes entanglement always sensitive to the size of the interval as 
we flow from the UV to the IR. This is strikingly different than in the absence of such anomaly. 
The dynamics of the probe that measures entanglement is described by the MPD equations rather 
than geodesics [26] , and our analysis shows that it exhibits the right features that explains why the 
system is not gapped in the IR. 

As highlighted in [26, 28], one of the interesting properties of the contribution of the gravitational 
anomaly to entanglement entropy is that the answer is now sensitive to the choice of Lorentz frame. 
To be more precise, if we rotate our normal vector by an S'0(1,1) transformation 

n“ —)> A^(r/')n^ (4.1) 
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with boost parameter r]', then, due to (3.53), entanglement entropy transforms as 


4e“’ ^ ^ A-)' ■ (“.2) 

This change is universal, and could be interpreted as one way to “measure” the anomaly from 
entanglement. This derivation holds both in the bulk and boundary theory, however, from the bulk 
point of view these manipulations are only consistent if there is a smooth connected worldline. For 
that reason it is crucial to explicitly construct connected solutions, which is a rather non-trivial 
task. 

We were able to construct only perturbative solutions to the MPD equations in two different 
regimes: small and large s. It would be very interesting to find exact solutions and hence to have 
a better understanding of S'gg as a function of s. Another interesting feature that we found 
was that any connected solution does not lie on a constant time slice, despite the fact that the 
background is stationary. This should be contrasted to the original Ryu-Takayanagi prescription. It 
is not clear if this feature will persist in the higher dimensional examples discussed in [27, 28, 29, 30], 
but it is worth further investigation. 

Along the lines of the construction in [36, 34], it would be interesting to build a monotonic 
function that controls the RG flow out of holographic entanglement entropy. From our explicit 
results in section 3.2.4, it is not evident how to build such quantity from However, 

a more interesting route is to modify the construction in [9]: for a QFT with a chiral anomaly, 
can we derive the constraints on cl^r along a RG flow via appropriate inequalities associated to 
entanglement entropy? We leave this question for future work. 
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is also the case that our bulk theories are not generically unitary, since TMG carries negative energy excitations. 
This feature makes it confusing whether we should even demand that any candidate c-function is monotonic. 
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A Conventions 


Throughout this work we will use the following coordinate system for the RG backgrounds: 

dz'^ rjijdx'^dx^ 
z‘^f{z) 

e 2 A ( u ) (^^^2 ^ , 

The relations between various definitions gives 

z = e“^ , , du = . 



In the limit IR limit, for 0 < n < 4, we have 


u ~ -A oo 

u ~ log(z) —)■ oo 
u ~ uo + > uq 


if 0 < n < 2 , 
if n = 2 , 
if 2 < n < 4 , 


and in the UV limit we have 


^^0, fiz) ~ 1 , u ^ z . 


B Linearized analysis 

In this appendix we derive the master equations for the metric and scalar fluctuations at 
level. The action is 

i[g, </>] = [R - 8{Vcj)f - v{<p)] + Igcs , 

and the equations of motion are 

Rfiu - \gf,uR + = 8df,(J)du(i) - Ag^yd°'(t)da(t> - \giJiuV{4>) 

^ j-l A 

and 

lev^,^ - = 0. 

To start we will keep the potential V{(f) arbitrary, and it will be specified as needed. 


(A.l) 

(A.2) 

(A.3) 

(A.4) 

linear 

(B.l) 

(B.2) 

(B.3) 
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We setup the linearized analysis by first defining 

9^^u = gjS +hf,u , 4> = 4>{o) + H ■ (B.4) 

Here {gjiJ, (j){o)) define the background solution, and {h^y,5(j)) are the small perturbations around 
the background. We use a radial gauge = 0. For the background we have 

[du^ + dwdw) , (B.5) 

and we will use null coordinates for the boundary directions; 

w = —t + X , w = t + X , d = dui , d = dyj . (B.6) 


When manipnlating the equations, the following identities are useful 

8((^(o))^ + -A{u) — (H(u))^ = 0 , + 2A + 8(^q = 0 

which are due to (B.2) at zeroth order, and the scalar equation (B.3) gives 



Here, and through out this appendix, prime denotes derivative with respect to u, and 

Fo = F((/)(o)) , 00 = 0n0o , A = duA{u) . 


(B.7) 


(B.8) 


(B.9) 


The linearized metric pertnbations are constructed as follows. Plugging in (B.4) in (B.2), we 
get that the (i, j) components of (B.2) are: 


1 

1 



2d‘^dh^w + dd‘^hu,w + + e^^dudh^w 

2dd'^hy,iu + + d'^dh^w ) + e^dudhyjw 


e^dudhww + 8e^0(o) d5(p = 0 , 
e^dudhww — 8e"^0(o) 064) = 0 ^(B.IO) 


and 


( ^6 (p^Q'jScj) H" 6 dy^hujuj 



(B.ll) 


The trace of (B.2) gives 

- - d'^hww + ‘2^ddhu,w + e~‘^^du [e^^duh^w) + 12(0(o) + A^(^q))64) + 40(o)5„(50 = 0 . (B.12) 
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However, this is not an independent equation: taking a derivative of (B.12), and after using heavily 
(B.lO)-(B.ll), one obtains 

- 1 

e~^du{e^duH) + duKu, = 0 (B.13) 

which is the equation obtained by linearizing (B.3). Hence the independent equations that describe 
the fluctuations are (B.lO)-(B.ll) and (B.13). Equation (B.12) serves as a constraint to fix inte¬ 
gration constants. In comparison with prior literature, our equations agree with [41], provided we 
set 5(p = 0 and the background metric is AdSa; note that ((9p)here = ~2(5p)there- 

B.l Einstein gravity plus scalar system 

Setting l//i to zero simplifies significantly the analysis. The linearized equations (B.IO) and (B.ll) 
simplify to 


e^dudhww = dKi , 
e^dudhww = dKi, 


8e^(j)^Q)54> -k e^duhww = Ki . (B.14) 

where Ki is an integration constant coming from integrating (B.ll). The components for 
and hiuu) of the metric are independent of dcj) and can be solved easily provided a warping factor 
e 2 A(y) jg specified. Using (B.14), the scalar equation (B.15) reduces to 


dli> + (a;2 - + IUeff(n)V' + 2Kie-^/2^(o) = 0 . (B.15) 


where V' = This is a Schrodinger equation (with a inhomogeneous term proportional to 

Ki) and effective potential 


lEeff 


—e 


„2A 


—e 


2A 


1 d^Vo 
16 

1 d^Vo 
16 


16 (^( 0 ))^ + e ^/‘^du{e^due ^/^) 


9 

4 


(A? + \a 


(B.16) 


^^Note that we use the term integration constant loosely as Ki is a function of x and t. We will use the notation 
where Ki is also expanded in Fourier modes in which case d and d can be viewed merely as numbers. We will use 
this convention for the remainder of the section and apologize if it creates any confusion. 
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In the IR, the effective potential behaves as 


VFeff ~ 7-^ ~ (B.17) 

[u - uor 

The homogeneous part of the solution for 6(j) is thus similar to that of a probe field in the IR. From 
this we conclude that the system is gapped if n > 2. 


B.2 Decoupling of the equations 


Here we briefly describe how to decouple the linearized equations for the metric perturbations: we 
will find a set of differential equations that only involve {hww, ^4’)- 

The key is to take two derivatives of the constraint (B.12) and then use repeatedly (3.10), (3.11) 
and (B.12) to eliminate {hisu,, h^w)- The resulting equation is 


[dp{e ^^5p) + e ‘^^Add - (/i^ - e 
+e“2^X" - e-^^A'Y” + 2ddX - 


—2A 

-^^H'2)] h'X - 

2A'e-‘^^ddY - - e-^^A''^) T' = 0 


(B.18) 


where 

X = 4e-2^(3())'('o)<5</. + 0'(o)5</>') , T = , 

and for compactness we introduced the radial variable p which is given by 


X' = dpX = -e^duX . 


In this notation the linearized Klein-Gordon equation (B.13) reads 


e-^^54>" + e-^^Add64) 


1 aVo 
16 902 


64> + 2e ^^"^0(0) h'^w — 0 • 


(B.19) 


(B.20) 


(B.21) 


Replacing (B.21) into (B.18), we get a fith order differential equation for <50. This equation 
is difficult to solve but one can study its IR limit. Taking this limit induces one simplification 
in the equation: the terms are subleading as <C the IR. The same mechanism 

responsible for the decoupling of the h^w hww is present here, the gravitational Chern-Simons term 
is dominating in the IR. However even upon taking an IR limit we cannot solve this equation for 
general n, but we can get an analytic solution for the special case n = 2. The solutions behave as 

(50~ (B.22) 

This behavior is essentially similar to that of the n = 2 case for a probe field [34]. It is thus likely 
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that the h^w and 6(j) helds have the same properties as a massless probe field and are gapped for 
n > 2. However, they only act as a source term for the diagonal components of the metric whose 
homogeneous equations still lead to ungapped modes. 


C U{1) gauge anomalies 


Another very interesting example where the boundary CFT can be made chiral and therefore stable 
against forming a mass gap is a bulk theory with the addition of gauge Chern-Simons terms. For 
simplicity, let us consider U{1) CS theories, contributing to the action the following terms 

Scs [ AAdA-^ [ AAdA , (C.l) 

Att J Att J 


where ki and kji are positive and A and A will lead to a left and right moving respectively U{1) 
current algebra at level ki, and kji at the boundary [51, 52, 53]. What is interesting about this 
scenario is that, as for the example with gravitational anomalies in section 3, the Chern-Simons 
action (C.l) does not affect the background (2.1). Therefore, in a holographic setup there is a 
potential ambiguity if the system is truly gapped or not. 

Before considering what is the effect of (C.l) in the bulk, we will review the CFT 2 reasoning 
behind the fact that the chiral anomaly is protected under RG. For any Lorentz invariant theory, 
ki — kji is conserved along the RG flow of any 1-1-1 dimensional theories that is a CFT in the UV 
limit. One version of this proof is as follows: current conservation dictates 

^ ^ (C.2) 


where the right hand side arises from quantum anomaly and T is a background gauge field that is 
coupled to the global symmetry. Due to this anomaly, the two point functions are 


{jwiw,w)j^niO)) = 


{ju,iw,w)jaiO)) = 


kLi\w\) 




a{\w\)'nww + b{\w\)ww 


re 


(C.3) 


kdH) 




where w = x + it and w = x — it. The above relations follow only from Lorentz invariance, which 
leaves at this stage the functions a(|ri;|), b{\w\), kL,R{\w\) arbitrary. Our goal in the following will 
be to restrict their dependence on |t(;|. 

We hrst note that for a CFT in the UV, approaches a constant as rc —)■ 0 and that 
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{jwjw) —^ 0. Note also that since and ju, commute when they are not coincident, it means that 


{jw{w,w)ju,{0)) = {ju,{w,w)j^{0)) . (C.4) 

Now we can use current conservation to relate the correlators given above. We are not coupling 
the global symmetry to any external field, so current is in fact conserved. Moreover, even if there 
is a non-trivial background field, the right hand side of the conservation equation (C.2) only leads 
to extra contact terms that would not contribute when the currents are not inserted at coincident 
points. Therefore we have 


{djuj{w,w)j^a{0)) = -{dju,{w,w)jw{0)) , 


(C.5) 


leading to 


dkLi\w\) _ ^^ a{\w\)r]^w + b{\w\)ww ^ 


(C.6) 


where we used (C.3) and (C.4). One can repeat the same exercise with replaced 

appropriately in (C.5) and obtain 


dkR{\w\) _ ^^ a{\w\ )r]^w + b{\w\)ww. 






(C.7) 


Now let us take w to be real, le t = 0, w = w = x. Then df{\w\) = 5/(|t(;|)|^g]R. Therefore, we can 
now take the difference between (C.6) and (C.7) to get 

[±{kL{x)-kii{x))]=0. (C.8) 

We have thus shown that k^ — kn is a conserved quantity as we increase x, equivalent to the effect 
of an RG flow towards the infrared. 


With an unbalanced kL — kn, we do not expect a gapped infrared theory. It is however clear that 
any solution in Einstein scalar theory continues to be an exact solution when the Chern-Simons 
gauge terms are included in the action. Moreover, in this case, the entanglement entropy formula 
is sensitive only to the metric and immediately we conclude that the log term is absent! 

The resolution of the apparent paradox is that the contribution of the C/(l) to the entanglement 
is invisible to the Ryu-Takayanagi holographic entanglement because the 1/(1)’s contribute to the 
central charge by an order one effect relative to the large classical contribution of the geometry; this 
effect which is subleading in the large N ^ ijG^ limit. To recover the effect of the current algebra, 
one has to take into account the quantum contribution to the path-integral of the Chern-Simons 
terms. 
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The quantum contribution of the Chern-Simons terms basically boils down to the contribution 
of the boundary modes. Let us consider the way it works in pure AdS space. The quantum 
contribution of the boundary modes to the action is given in terms of modular functions. 

We are in particular interested in the quantum contribution of the boundary modes to the 
entanglement entropy. To deal with the replicated geometry, we could play the same trick as in 
[54] and make use of the hyperbolic slicing. The replica index n becomes the temperature of a bulk 
hyperbolic black hole, and controls the periodicity of the time coordinate. One can compute the 
boundary partition function using the Cardy formula, where the time direction has periodicity 

^ = I3n = 27rn , (C.9) 

and spatial volume is infinite, but has a cutoff that takes the form 

L = 21og(|), (C.IO) 

where R is the size of the interval, and 5 a short distance cutoff from the boundary of the inte¬ 
gral.Together, this gives 


log Zn = log(trexp(-,diL)) = 
and finally the entanglement entropy is given by 


7r(cL cr) L 


12 


/3 ’ 


See = 


n Tn log Zn - Ti log Zi 


1 — n 


Ti 


CL + Cr .R, 


(C.ll) 


(C.12) 


For each copy of the gauge field with positive (negative) level, it contributes to a left (right) 
moving chiral bosonic mode at the boundary. This contribution to the entanglement entropy is 
thus clearly subleading in 1/iV^ compared to the gravitational sector. 

The point is that since the Chern-Simons terms are topological, it is insensitive to the internal 
geometry so long as it is regular (in this case it is important to regulate the curvature divergences 
mentioned in section 2), and so the only contribution to the quantum action comes from the 
boundary, which stays the same, unaffected by the holographic RG flow induced by the neutral 
scalar field U{1). It contributes to a log term confirming that the infrared is not completely gapped, 
although these massless modes are not visible in the leading large N limit. 
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